Primordial non-Gaussianity estimator: the inhomogeneous noise effect 
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Since the inhomogeneous instrument noise can produce extra non-Gaussianity in the CMB 
anisotropy, its effect should be carefully subtracted in the primordial non-Gaussianity estimation. 
| We calculate the probability distribution function of the CMB anisotropy for local type of non- 

, Gaussianity, from which the optimal estimator in the general case (inhomogeneous noise and cut 

■ sky) is obtained. The new estimator obtained here is different from the popular one, since the 

' inhomogeneous noise and cut sky effects are completely accounted. The CMB anisotropy in the 

new estimator is noise weighted. The noise weight is different from that used by WMAP Group in 
their 5-year data analysis. Although it is still difficult to calculate the new estimator rigorously, for 
the case of the slightly inhomogeneous noise, there exists a series expansion method to compute the 
new estimator. Each order in the series is suppressed by two factors, (^ — 1) and ^tfst, which make 
the method feasible. Through the Edgeworth expansion we can generalize our discussion to other 
types of non-Gaussianity. 

6 
& 

I. INTRODUCTION 

The standard inflation predicts a flat universe with scale invariant and gaussian random primordial fluctuation, 
£> ' which is confirmed by the Cosmic Microwave Background (CMB) observation [lj. Despite its extreme success in 
modern cosmology, there are too many inflation models which all give the above predictions. However, with increasing 
data, it is possible to distinguish those models by some subtle features, for example, the primordial gravitational waves, 
the running power spectrum, the isocurvature perturbation and the primordial non-Gaussianity [2]. Among them the 
primordial non-Gaussianity is an important one. 

The primordial non-Gaussianity can be quantified by the 3-point correlation function of curvature perturbation 
<£>(fc) (equivalently by the bispectrum B(ki, ki, k 3 )), 

o 



^1- 
o 



PACS numbers: 98.70.Vc, 98.80.-k 



< $(fc!)$(fc 2 )$(fc 3 ) >= (2n) i S J (k 1 + k 2 + k 3 )B(k!, k 2 , k 3 ). (1) 



k> I Acoording to [2| the primordial non-Gaussianity can be roughly divided into two types: local type and equilateral 
type, "local type" means B(ki, fa, k 3 ) is large when one of the momenta is small compared with the other two, while 
"equilateral" type means B{k\, ki, k 3 ) is large when the three momenta are of the same order. In the real space the 
local type of non-Gaussianity can be represented as: 

$(x) = $ L (x) - f NL ($l(x)- < $!(£) >), (2) 

where $j, is a gaussian random curvature perturbation and /jvl describes the strength of the primordial non- 
Gaussianity. 

The standard inflation model predicts only a small non-Gassianity 0, which is beyond the experimental limit. 
However, this is changed in multiple-field models, which is investigated in @, 0] and many other papers. Besides 
this, non-local inflations and ghost inflation theories can also produce large non-Gaussianity [a. M, Ho|. Even in the 
single field inflation, it is possible to generate large non-Gaussianity [Til fl2||. 

So the non-Gaussianity measurement is important to constrain the inflation dynamics. The general properties of 
primordial non-Gaussianity in the bispectrum measurement are discussed in However, the direct measurement 
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of the bispectrum is very time-consuming for WMAP and future data. A fast local type of non-Gaussianity estimator 
is proposed in [l4[ and it was applied to the 1-year WMAP data to constrain the non-Gaussianity in [ia |. The fast 
estimator proposed in [3] is originally constructed on Wiener-filtered map s. Later on, [l6[ proved that it is optimal 
in weak non-Gaussianity and full sky case. But it is already realized in [14| that the inhomogeneous noise and cut sky 
play an important role in the non-Gaussianity estimation, since from them extra non-Gaussianity can arise. In [171 ] 
a linear term is introduced into the estimator to correct this inhomogeneous noise and cut sky effects. The complete 
estimator can be written as: 

Qin(a) = — [(< a-l 1 rn i ai 2 m 2 ai 3 m 3 >f NL =1 )C( a )i imi ,i 4 m 4 C( a )Z 2 m 2 ,i 5 m 5 ^ (a)l 3 m 3 ,l 6 m 6 a hm 4 ai 5 m 5 ai 6 m 6 

— 3(< ai imi ai 2m2 a l3m3 > f N L =l)^ ( a )i imi: i 2m2 ^ (^hmsJimi 01 ^™-^^ (3) 

where C7 a ); imi! ; 2m . 2 =< a/ 11Tll a; 2m2 > and S is the normalization factor. In real calculations, it is difficult to obtain 
CTn if inhomogeneous noise and cut sky are considered, since C( a ) is a huge matrix when I is large, it is impossible 
to find its inverse directly. But once C~^a can be quickly calculated, a feasible algorithm based on this estimator 
can be executed [l8|]. Nowadays, is approximated by ^f 2 - fLU [T7L Il9j. which is easy to calculate. After this 

approximation the estimator becomes, 

/ \ _ }_ \ ^ r (< ai 1 m 1 0-l 2 m2 a l3m3 >f NL = l) 
£Un\a) — ~ I C C C a hm 1 a,l 2 m 2 ai 3 m 3 

(l,m) h h h 

Q (< Qj 1 m 1 Q; 2 m 2 Qi 3 m 3 > f N L = l) n 1 (A\ 

n n r< (a)2i"ii,i 2 «i2 a '3"i 3 J- 

Besides in the non-Gaussianity estimation, the same problem also arise in power spectrum estimation [20j and CMB 
lensing detection plj . 

Although the inhomogeneous noise and cut sky effect is already considered in the estimator J3J 2]), it is incomplete, 
since the noise and sky effect enter into the estimator only through the spherical harmonical coefficient et/ m . The 
aim of the present paper is to give a complete treatment of the inhomogeneous noise. We calculate the probability 
distribution function (PDF) of CMB anisotropy in the inhomogeneous noise and cut sky case. From the PDF a new 
estimator is obtained. For slightly inhomogeneous noise it is possible to compute the estimator order by order. Each 
order is suppressed by two factors: — 1) and ^=&, so the series converges. 

This paper is organized as follows. In section II, the PDF for the CMB anisotropy is considered, from which we 
obtaine the optimal estimator for local type of non-Gaussianity. In section III, we expand the new estimator, and 
get a series expansion method to deal with the inhomogeneous noise. Through a detailed analysis we show that in 
what condition our estimator reduces to the popularly employed one and how our method reflects the inhomogeneous 
effect. In the last section, we first generalize the estimator to other types of non-Gaussianity, then we also discuss 
the slightly inhomogeneous condition and the noise weight. Appendix A is devoted to the relations appearing in the 
radiative transfer process, while Appendix B is devoted to the normalization factor. 

II. THE BISPECTRUM ESTIMATOR 

The observed CMB anisotropy is composed of several components: the primordial temperature fluctuation, the 
foreground emission, secondary effects and the instrument noise. In this paper we do not consider the foreground 
emission and secondary effects, since the foreground emission has already been subtracted [22j and secondary effects 



such as point sources and CMB lensing are carefully studied in [23], |2J, |25( . Therefore each pixelied CMB data can 
be regarded as superposition of the CMB signal and the instrument noise, 

5T, = 5T? + 5T>\ (5) 

For the ith pixel, the total PDF is just a simple convolution of the PDF of the signal and the PDF of the noise [26j |: 

f(STi) = [ dST?dSTpf,(STf)f n (STp)S D (m - 6Tf ~ 5T?), (6) 



where f s is the PDF of the signal, f n are the PDF of the noise and Sj) is the Kronecker delta function. The instrument 
noise for each pixel is always assumed to follow the gaussian distribution with a variance <ii , 

f (ST n ) 2 

/J«r) = -=exp[-i-i^], (7) 

V /27T(7/ Mi 
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where <7j varies from pixel to pixel. For a given primordial curvature perturbations ^(a;), the PDF of the signal 
f s (ST?) can be written as: 

f s (6T?) = f d N <P5 D (ST° - f r 2 drJ2^(r)biJ2^(r)Y lm (n l ))f(^), (8) 

•* ^ l m 

where 6/ is the window function which reflects the beam and pixel smearing, /($) is the PDF of the curvature pertur- 
bation and ai (r) is defined by EQ (|A4|) . (r) is the spherical harmonic coefficient of the curvature perturbations at 
a given distance r and n, is the unit vector pointing to the %th pixel. For the local type of primordial non-Gaussianity 
(J2J, the PDF /($) is given by: 

/(*) = / ^$ L ^(d>(f) - * L (f) - /^[«i(50- < >]) C ^^^ L] ' (9) 

where the functional /(<&) is written in compact vector notation. In the spherical harmonical space, C$ is given by 
EQdXTJ. Substituting EQ 0, ©, and © into EQ ©, we obtain the one point PDF for the i th pixel, 

S{m) = f d N ^d N <i> L dST l n S D (5T l ~ I r 2 dr^ai(r)biJ2®lm(r)Yi m (ni)-5T?) 

l m 

' : exph^!]^($(f) - ^(af) - (*)" < >])^4™==^- (10) 



The above formula can be easily generalized to the united PDF for more pixels, 

f(ST u --- ,ST m ) = J d N $d N $ L Y[{d5T- l 8 D (STi - J r 2 dr ^ a/(r)& ; ^ <J> ;m (r)y Jm (ri. ( ) - ST") 



i=l " I 



Today's CMB constraint indicates that the primordial non-Gaussianity is weak, so we can calculate the PDF ([9]) to 
C(/jvLA t ) (a* =< & 2 (x) >) and safely ignore higher order terms. This is already done in [l6(. Here we just present 
the final result, 

/($) = exp[-i^ J rldnr^^MD^in^M 

l,m 

+ /nlJ2 /^r 2 2 dr 2 £^^ (12) 

(l,m) J 

where D~ 1 (r 1 ,r 2 ) is defined by EQ (T£3j) . and g^' 2 ^ 3 2m is Gaunt Integral defined by EQ (|B2|) . Then substitute this 
result into EQ (fTTj) and do the ST integral, the PDF (JTTJ) is simplified to a Gaussian functional integral, 



with 



/d d d 
rXdr^ld^Q 1 ^^— — D^ 1 (n , r 2 )— — - — — -] 

{ [ d »* lmd »$ tm (l[—^)eM-j:S] 

exp[-$*C$ + J*$ + J$* + A*$ + (13) 



2 '2 

C irori(Iror) ' = r -^{DT\r,r')5 ul 5 mml + Q Z (r)b,Q z - (r N^ [lm)l ) (14) 



lm,(lm) 



E 



-M(ni) 



(15) 
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A lm (r) = r ^&l STY^u^M^), (16) 

where M(rii) is the mask function, which take value unity when i = l,...,m and zero in other case. Here and 
in the following we use the following abbreviations, $*C$ = J drdr E(imi'm') ®im( r ) C lmr,(lmry ®i' m ' = 
I drJ2im J*m( r )®lm( r ) an d = / drJ2im ^fmM^mM- After performing the integral over <!> and dropping the 
irrelevant constant factors, we obtain 

\gf = A*C- 1 A + f NL ]T / r^r^r^^ (17) 

(l,m) J 

It is obvious that performing the Gaussian integral is equivalent to simply substituting $; m by C~ 1 A. Acoording to 
16], the optimal estimator of /jvl can be taken as: 

< T ) = | E / ^^^^(C- 1 ^ (18) 

(Z,m) 

where S is the normalization factor. 

EQ (|18p is our first main result. In deriving formula (|18[) we do not require the condition of full sky and homogeneous 
noise, so EQ (| 18|) can be applied to cut sky and inhomogeneous noise. It should be noted that even if one only considers 
the cubic term (we do not consider linear term for the moment), EQ p8[) is also different from EQ ([3]). A in EQ 
(fTHj) corresponds to the ordinary spherical harmonical coefficient ai m in EQ ([3]), but it is noise and sky weighted. C 
in EQ (fTgj) corresponds to C( Q ) in EQ (J3J), but C is a function of r while C( a ) is not. This is because only the PDF 
of curvature perturbation is approximately considered to /nlH order, the cut sky and inhomogeneous noise effects 
for each pixel are rigorously accounted, while in writing down EQ the cut sky and noise effects are only partially 
accounted through the ordinary spherical harmonical coefficient ai m . 



III. THE INHOMOGENEOUS NOISE 



Although EQ (fT8|) is a general estimator in the case of inhomogeneous noise and cut sky, it is difficult to calculate the 
inverse C _1 A, which becomes the main problem in the real computation [l8j . Luckily there exists a series expansion 
method to deal with the inhomogeneous noise when the noise is slightly inhomogeneous. In order to establish the 
conventions, we first consider the case of full sky and homogenious noise and then turn into the case of cut sky and 
inhomogeneous noise. It must be stressed that the method presented in the following cannot be used to deal with the 
cut sky effect even there is no noise. In the subsequent section the window function 6/ has been absorbed into the 
definition of a;,/3;,C;, therefore when a>i,/3i,Ci appear, they mean aibi,(3ibi,Cibf. 



A. Full sky and homogeneous noise 

In the case of full sky and homogeneous noise, the noise level is the same for each pixel (denoted by <r), so the 
estimator (jT8)) will be greatly simplified. First note that in this case A -1 , C, A reduce to 

r-l 



N lm,(lm)' ~ 2 4tt V ^mm' > ( 19 ) 

N pix 



r cm (t) 

A lm (r) = 2 l \J a lm , (21) 

1(7 m~ 



where a; m = X^=:T ^Zm(^)<^ ^ s the ordinary spherical harmonical coefficient. It is easy to find the inverse C 1 



(C-i) lmr = 2(A(r,r') - A(r ff (22) 
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where C\ ot — Ci + <^ 2 ^- is the total angular power spectrum. Then C 1 A becomes 



(C 1 A)lmr = ^Qtaira, (23) 

which gives the optimal estimator, 

£ / Tl _lv /" Ja^iWs a h(r) „ ft»W . 4W n , 9zn 

n \ « Zi 



0, 



This is the estimator used in [14 ] . So for full sky and homogeneous noise, our estimator (fT5)) reduces to the ordinary 
one, as is expected. 

B. Cut sky and inhomogeneous noise 

In the case of cut sky and inhomogeneous noise, the weight (including the sky mask and the noise level) is different 
for each pixel, so it is impossible to calculate the inverse of C directly. However in the case of slightly inhomogeneous 
noise, there exists a simple way to calculate C _1 A. In this case, we split the matrix TV -1 into two parts: 



i=\ 

1 

j-2 ^ °W °mm +iV i m: (i m )' 



Su'S^+N- 1 * (25) 



with 



1 1 N 



a 2 N^a, 

2—1 1 



- E KifcWmW&Mfc) - 1). ( 27 ) 



N- 1S 



Zm, (Zm) 



i—i 



where TV is the total number of the unmasked pixels. The first term in the last line of EQ (|25|) represents an average 
effect, and the second term is a small perturbation about the average. Correspondingly, A can also be split into two 
parts: 

C - r —(D-Hr A I ^Mll vr s r*ai(r)r 2 ar(r') 1S 

^lmr,{lmr)' - r, \ u l \ r i r ) T 2 4tt WW+ 9 JV im,«m)' 

— ° )lmr,(lmr)' + (^)imr,(lmr)' ■ (^8) 

The first term is just EQ (|2TJ)) . its inverse is given by EQ ([2"!?]) . Compared with the first term, the last term is small, 
so we can calculate the inverse C -1 order by order, 

c -l = C {0)-1 + C {1) + C {2) + . . . + C (n) + . . . j ( 29 ) 

where C" is given by = (-^(C^^B^C^- 1 . From EQ (28]), C^^B is represented as 

(C^B) lmrSmr)> = ^^ y r\(r), (30) 

and C^ -1 ^ is given by 

(C (0) -^W = ^(* 2 ^&i m ), (31) 
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with 

hm = J^Y^iH^MiHi). (32) 

i=l » 

Thus, combining with EQ (|29]). EQ (J30J) and EQ pT]). (C _1 A) iml . can be calculated by the scries: 

(C^A)^ = (C^A)^ + ■■■ + (-l) n ((C^- 1 B) n C^ 1 A) lmr +■■■ 

_ A(r), 2 4tt _ #(r) ^ 2 4tt 1S C ix 2 4tt 

— /^tot v" jvr --ftot W ^Im^hmi) r tot \" at °hmi) 

1^1 J-^pix L'; j m J-^pix ii Jv pix 

AO) ,_a 4tt Ar _ lg , C fl fJZ 4tt Ar _ lg A / jJli, n , 

^ /^fo* ^ AT ly lm,l imi ) x CT A T JV Zimi,Z 2 m2'' r'tot K a l\r y i2iri 2 ; + "- 

1 l\m\l2m 2 61 12 

/3;(r) 2 47T ■r-^ 

= T^T " (Oim + } ^lm,l'm' b l'm') 

I m 

AM^ m x 

In last second line we formally define the matrix F which will be used in Appendix B and in the last line we define 
a quantity T; m which is analogous to the spherical harmonical coefficient a/ m . With Tj m the estimator (|T8l) can be 
written in a compact form, 

e (T) = iv [r 2 drG 1 ^ ^1 T , M1 T , ^-T, (34) 

(l,m) ll 12 13 

This is our second main result. Compared with the original one (|18p . now the estimator has the same form as the 
ordinary one (|24)l . the only change is substituting <z; TO by T; m , which is calculated by the seires (|33|) . Each order in 
the series can be computed as easily as a; m . In the following two subsections it will be shown how Ti m takes the 
inhomogeneous noise effect into account. 



C. Full sky and inhomogeneous noise 



In this subsection we consider the case of slightly inhomogeneous noise but full sky. Under this condition, M(n;) 
in N^ (lm)l (H7J), k m ([32D all equal unity and a ^ is determined by 4* = — |- J^fj-i ^- From EQ @, we see 



_1 Q 

that there are two factors to suppress the series. The first one is the factor ^ — 1 appearing in N^^., , and the 
other is the factor -S^. In the large scale region the latter factor is close to unity, however, with I increasing, where 

noise begins to dominate, it becomes smaller and smaller. The first factor does not vary with I, it is determined 
by experiment and small as long as the slightly inhomogeneous condition is satisfied. Therefore, the effect of the 
inhomogeneous noise can be accounted order by order through (I3"3"f . 

In order to show what T\ m is compared with o; m , we keep C~ 1 A (|33[) to first order 

(r -i A \ _ A(0 T 

AO) _2 4 *" /> V- AT -is , c i' w 2 47r , « /ot-x 

I m ' 



and split 6; m into two terms 



&im = 2 47r (aim + dim), (36) 



with a ; .,„ = j^-J^f^r Y hn(™i) ST i and rf /m = w^Ya^i ^ m (^)^(fr - !)• Compared with a im , d im is a small 
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quantity. Pluging ai m and di m into EQ (|35() and ignoring di m in the last term in EQ (|35[) . we obtain 

4tt JV " i '" cr 2 4tt JV " iB cr 2 C/ 

T im = a lm + —— E Ylm(ni) ST i(-2 - !) + T7 - XI ^mfaX^ ~ *) X y /'m'(^)7^t a i'm' 

iv p" <=1 iv p" i=l l> m > u l' 

= Aim + "T7 X] y imK)(rS ~ 1 ) X] y ''"' ^ rtot a l'm' • ( 37 ) 

It is easy to see that the dominant part of TJ m is the ordinary spherical harmonical coefficient a/ m . Besides this term, 

jv — — 2 

there is another small term which is suppressed by factors c t„r and 2j — 1. This term reflects the inhomogeneous 

cr 2 4,r 

noise. It should be pointed out that the factor c tir appears in the difference T; m — a; m , while the factor 

appears in the series (|33|) employed to calculate T/ m . So in small I region, where the CMB signal dominates, T/ m is 
close to aim, but the series ((33)) converges slowly. While in the large I region, Ti m is different from a; m , but the series 
converges rapidly. 



D. Cut sky and homogeneous noise 

In the previous section the case of full sky but inhomogeneous noise has already been discussed. Here we mainly 
discuss the case of cut sky but homogeneous noise. In this case, the second term of the series (|3"3")) becomes 

4?r AT -1S >, C h / 2 4?r ■ \ 
(imi ii v 

= tt— X y 4(^)(^K) - 1) E *W*)r5t aT" X ^(n.O^-M^), (38) 

1* pix - n , ^7. ^ pix ■ -, 

where (o' 2 ^ ZI -6; imi ) = X)j=T ^* 1 m 1 (ftj)STjM(fij) represents the temperature fluctuation of the unmasked region. 
In the small I region, where the signal dominates, -M^r rs 1, EQ ([38]) becomes 

4ir Np >* 4tt N ' ix 4tt 

' E y 4(n0(^(ni) - 1)t^ E 5T J M{n J )5 Dv {-—r 1 



47T ^ 

E ^4(^)<mM(n i )(M(n i ) - 



/v. 



ij=l 



= (39) 

In the large I region where the noise dominates, is small, so EQ (1551) is greatly suppressed. Therefore, in this 

i 

case the second term in (1331) is also small compared with the first term. 

However, we must stress again that the series ([33)) can only be employed to treat the effect of noise, but failed to 
account for the effect of cut sky. This is because in the zero noise limit, -M^ strictly equals unity for all I, all terms 

in the series ((33)) are zero except the first one. 



E. The normalization factor 



The last problem is the normalization factor S. S in EQ ([34]) can be calculated formally order by order. This 
calculation is presented in Appendix B, but the form presented there is not suitable for numerical computation. So 
we need a fast algorithm to compute the normalization factor S. However the advantage of estimator ([34)) is that the 
normalization factor needs to be computed only one time, while the part left which needs to be calculated repeatedly 
(to obtain the confident region of /nl) is easy to compute. 
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IV. DISCUSSION AND CONCLUSION 

(1) Until now we only consider the estimator for local type of non-Gaussianity. It is straightforward to generalize 
our results to other type of non-Gaussianity. Generally speaking, the PDF of non-Gaussianity curvature perturbation 
can be expressed by the Edgeworth expansion. The Edgeworth expansion of a multivariate PDF can be found in [13] 
and [28| . EQ (7) of [27j indicates that the PDF of the curvature perturbation can be represented by, 

/($) = (1 - X)/ ($), (40) 

where /o( < i ) ) is the Gaussian part of the PDF and X is a small perturbation which is function of (the exact form 
of X for general non-Gaussianity can be found in [27j). Following the same steps of section II, where the key point is 
to substitute $>i m in X by C~ X A, we obtain the estimator for general non-Gaussianity, 

e (T) = ±X(C~ l A), (41) 

and the linear term is already naturally included in the estimator. 

(2) Now let us discuss the slightly inhomogeneous condition. Here we simply take the WMAP 3-year V2 band as an 

2 2 

example. The noise level for each pixel is determined by the observation times af — N , 2 — j^° Jf , where <7q is 
the rms noise per observation. From EQ (|26|) we see that the effective observation times corresponding to the average 
noise level is N e ff = -k J2i=i N bs(i)- Since there are about 4% pixels with N i, s > 1200, 1.8% pixels with N b s < 400 
and for most of pixels the observation times is distributed in the region [500, 1000], the factor (^ — 1) = ( A /v" / ^' > — 1) 
will be small for most pixels. Those pixels with too large or too small observation times can be masked out in the 
non-Gaussianity estimation. 

As a matter of fact, there is no requirment that o 1 must take the form of (|26|) . It can be any values. However in 
order to make the series calculation efficient, EQ (|2T)|) is a good choice. There are two factors to affect the choice of 
a 2 , (2tt — 1) and If <J 2 is too large or too small, the series (|33|) will converge slowly. 

(3) In preparing our paper, WMAP Group release their five-year data 29]. In order to deal with the inhomogeneous 
noise, they adopt the average noise and "combination signal-plus-noise" weight explicitly in EQ (A6), EQ (A18), 
EQ (A19) and EQ (A27) of [2^]. Using the symbol conventions in our paper, they are a\ = X^=:T °f -^(^)j 
(^+^) 3 = ^EiLT(^g^) 3 . Wiru) = and a lm = ^ E^T ^n^M^^^t (we already 
reexpress their equations for the present purpose), where cr 2 mb = 4- X)z(-^+l)Ci (remember that 6/ is already absorbed 
in the definition of Ci). In \2Q\ a\ is used to calculate the total angular power spectrum C\ ot = Ci + j^-af, while 
<72 is the average noise appearing in the "combination signal-plus-noise" weight, they are not the same. Besides this, 
the noise weight is also different from the way used in our paper. We do not know why they use such a "combination 
signal-plus- noise" weight, but our discussion indicates that the inhomogeneous noise effect will be accounted better 
by the method presented in previous section, where the inhomogeneous noise is treated seriously. 

In conclusion we have studied the inhomogeneous noise effect in the CMB non-Gaussianity estimation. First, we 
calculate the CMB anisotropy PDF for local type of non-Gaussianity which gives an optimal bispectrum estimator 
(|18[) . The estimator obtained in this paper is different from the ordinary one, since the cut sky and inhomogeneous 
noise factors are completely accounted in the PDF calculation, which provides a framework to study the effects of 
the cut sky and inhomogeneous noise in the CMB non-Gaussianity estimation. Then, we propose a series expansion 
method to calculate the new estimator. The final result is just to replace ai m which appears in the popular estimator 
(|24|) by Ti m . Tim can be calculated by the series (|33f as easily as ai m - To the zeroth order, Ti m reduces to ai m and 
the estimator ([3~4")l becomes the usual one. Last, we discuss how to generalize the estimator (|18p to other type of 
non-Gaussianity, the slightly inhomogeneous condition and the noise weight. 
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APPENDIX A: RADIATIVE TRANSFER 

This appendix summarizes the relations which appears in the radiative transfer process. The contents are mainly 
based on section IV of 

In the spherical harmonic space, a given curvature perturbation <f>(x) is represented by $; m (r) . The covariance 
matrix of $z m (r) can be calculated by 



{®*i imi { r i)®i2m 2 {r2)) = 5 hh 5 mim2 D h (r 1 ,r 2 ) (Al) 
= S h i 2 S mim2 - I k 2 dkP(k)ji 1 (kr 1 )j h (kr 2 ). (A2) 



(A3) 



7T _ 

where P(k) is the power spectrum. Symbolically one defines the inverse of Di(r\, r 2 ) as 



Introduce the symbols: 



ai{r) = - / fc 2 dfcjKMAz(fc), (A4) 



7T 



AW = \ J fc 2 dfcP(fc)^(fcr)Ai(fc), (A5) 
where Aj(fc) is the radiation transfer fcuntion [301 ] . Then we obtain the useful formulas: 

fr 2 dra l (r)P l (r) = C l . (A6) 

and 

«i(ri) = J rldr 2 Di l {r u r 2 )f3i{r2), (A7) 

A(ri) = J rldr 2 Di{ ri ,r 2 )ai{r 2 ). (A8) 

Finally we define the reduced bispectrum for local type of non-Gaussianity, 

Khh = 2 / r 2 dr[a h (r)A 2 (r)j8„ (r) + & (r)a Ja (r)$ 3 (r) + h {r)/3 h (r)a h (r)]. (A9) 

When the window function is considered, EQ (|A6[) . (|A8[) . (|A8|) . (|A9|) do not change except the replacement a/(r) — * 
ai(r)bi, f3i(r) -> /9/(r)6j, C; -> C;6f , 6; li2 ; 3 -> b h i 2 i 3 b h bi 2 bi 3 . 

APPENDIX B: THE NORMALIZATION FACTOR 

In this appendix we give the normalization factor formally. First, the 3-point correlator of <5Xi can be calculated by 

< STiSTjSTk > = < (dT/ + 6T?)(6T? + ST?)(6T% + 5T£) > 
= < 6T?5T?6T% > 

> 6z 1 ^ 2 6i3yi 1 m 1 (ni)5 / i 2 m 2 (nj)Fi3 m3 (nfe) 

= &m 1 m 2 m 3 ^li2i3^rm 1 (^)^ 2 m 2 K)^ 3 m3("fc) I ( B1 ) 

where bi 1 i 2 i 3 is the reduced bispectrum defined by EQ (|A9|) (the window function is already absorbed in the definition), 

O l rk[t 2m3 = I d 2 nY lumi {n)Y l2 , m2 {n)Y htm M- (B2) 



3 

Gm[ 2 m 2 rn 3 is tne Gaunt Integral 
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When deriving EQ (|B1[) , we have used the follwing formula 



i 3 

Then the 3-point correlator of bi m can be formally represented by 



bhhh- (B3) 



< b hmi b l2m2 h 3m3 > = E Yl * mi } nt) F ^>^ *W"0 M{ft . )M{Hj)M{fik) < gTiSTjST). > 



ijk^l CT * °i a k 

/-^i mi m 2 m 3 -^--^ cr? ^ 

I'm' *=1 2 J =1 J 

^ ^; m3 (n fc )y ;Xi (n fc ) 

2^ 2 JL - 2 — M ( n fe) 

fe=i °* 

= ij' 1 ' 2 ' 3 fB41 



mim2m3 ! 

With the 3-point correlator of bi m , the 3-point T; m can be formally written as 
< T, Ti Ti > = ff' 1 ' 2 ' 3 

rl'hl 



t m 

+ [Vf, / F, ,/ ,i^ 3 , + cycZ.l 

L Z_^ IimiA" 1 ! l2rn 2 ,l 2 m 2 m \m 2 m 3 y J 



4- \ F 1 1 F 1 ' F 11 H (B5) 

/ / hm 1 ,l 1 m 1 l 2 m 2 ,l 2 m 2 l 3 m 3 ,l 3 m 3 m'^m'^m',' ^ ' 

I rn 



Where the matrix F is defined in EQ (|33|) and H is defined in EQ (|B4j) . So the normalization factor S can be 
formally represented by 



9 - V / r 2 drG 1 ^ ah(r) Pi 2 {r) (3 h {r) 

°~ / ' aly m 1 m, 2 m 3 Qtot (Jtot Qtot ^ J -hmi- L l 2 m2- L hm 3 ^ ■ \ D <J) 

(l,m) ll 12 ' 3 
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